Frustrated pentagonal Cairo lattice in the non-collinear antiferromagnet Bi 4 Fe 5 Oi3F 
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The crystal and magnetic structures and underlying magnetic interactions of Bi4FesOi3F, a model 
system for studying the physics of the Cairo pentagonal spin lattice, are investigated by transmission 
electron microscopy, low-temperature synchrotron x-ray and neutron powder diffraction, thermo- 
dynamic measurements, and density functional band-structure calculations. The crystal structure 
of Bi4Fe50i3F contains infinite rutile-like chains of edge-sharing FeOe octahedra interconnected by 
the Fe2C>7 groups of two corner-sharing FeCU tetrahedra. The cavities between the chains are filled 
with the fluorine-centered Bi4F tetrahedra. The Fe 3+ cations form pentagonal units that give rise 
to an unusual topology of frustrated exchange couplings and underlie a sequence of the magnetic 
transitions at Ti = 62 K, T 2 = 71 K, and T N = 178 K. Below Ti, Bi 4 Fe 5 Oi 3 F forms a fully ordered 
non-collinear antiferromagnetic structure, whereas the magnetic state between Ti and Tat may be 
partially disordered according to the sizable increase in the magnetic entropy at Ti and T 2 . There- 
fore, Bi4FesOi3F shows the evidence of intricate magnetic transitions that were never anticipated 
for the pentagonal Cairo spin lattice. Additionally, it manifests a sillimanite (AbjSiC^-based ho- 
mologous series of compounds that feature the pentagonal magnetic lattice spaced by a variable 
number of octahedral units along the rutile-type chains. 

PACS numbers: 75.25-j, 75.50. Ee, 61.66. Fn, 75.30.Et 



I. INTRODUCTION 

Magnetically frustrated systems entail a pattern of in- 
teractions that can not be satisfied simultaneously.^ 
These patterns typically involve equivalent antiferromag- 
netic (AFM) interactions arranged in closed loops with 
an odd number of bonds, as in a pentagon or in a trian- 
gle. As one moves around the loop, the spin flips on each 
bond, and the translation symmetry is violated because 
the full turn around the loop leads to the inverted spin 
direction. Therefore, in frustrated systems spins can not 
develop a simple collinear AFM order, instead more com- 
plex magnetic ground states appearP Finding appropri- 
ate model systems for experimental studies of frustrated 
magnetism is a long-standing challengeP^D The availabil- 
ity of model systems is crucial for addressing theoretical 
problems and raising new questions to theoretical as well 
as experimental research. 

The loops with an odd number of bonds are typically 
triangular. P lanar triangular spin lattices have been 
widely studiecP^IMl m a t attempt to achieve a spin- liquid 
state that lacks long-range magnetic order down to zero 
temperature. When triangles form corner-sharing tetra- 
hedra instead of tiling the plane, the system enters a 
so-called spin-ice state that facilit ates experimental ac- 
cess to magnetic monopoles! 1 * 2 * 6 * 9 ! Thanks to a different 
topology, spin lattices based on pentagons could lead to 
a range of hitherto unexplored properties and ground 
states. However, experimental and even theoretical con- 
siderations of such lattices remain challenging because 
the five-fold symmetry of an ideal pentagon excludes the 



periodicity of the lattice. 

Pentagonal units do not form a regular tiling in the 
plane. To build a tiling, the pentagons have to be 
distorted and should, unlike the triangles, involve non- 
equivalent bonds. One of the least-perimeter tilings al- 
lowed for pentagons is known as the Cairo lattice. It 
entails the coupling J\ between three-vertex sites, and 
the coupling J2 running from each three-vertex site to 
the surrounding four- vertex sites (Fig. [lj. Recent theo- 
retical work puts forward a number of unconventional 
ground states, including spin nematics, an d oth er in- 
teresting properties of the Cairo spin latticeP^^ These 
predictions remain highly challenging for experimental 
verification. In crystals, the pentagonal arrangement of 
magnetic ions is very rare, which is partly related to the 
fact that periodic crystals can not have five-fold sym- 
metry. Presently, the only known material prototype of 
the Cairo lattice is L^Fe^g having mullite-type crystal 
structureP^HMI 

While theory readi ly sp ans the broad range of the 
J 1 and J 2 parameters) 10 1 11 1 Bi 2 Fe40g represents a single 
point on the effective J\-J% diagram.^ Experimental ac- 
cess to the diverse magnetism of the Cairo lattice requires 
a range of compounds with different J\ and J2. In the 
present paper, we demonstrate that a deliberate modi- 
fication of the sillimanite-type (A^SiOs) structure may 
lead to interesting materials showing exotic magnetism of 
the Cairo spin lattice. Specifically, we report the crystal 
structure and magnetic behavior of the hitherto unknown 
oxyfluoride B^FesOisF that undergoes a sequence of in- 
tricate magnetic transitions and reveals unconventional 
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FIG. 1. (Color online) Cairo pentagonal lattice (left panel) 
and its derivative version in Bi4FesOi3F (right panel). Note 
that the ideal lattice comprises only two couplings, Ji and 
J2, whereas in Bi4FesOi3F J2 splits into J a ti and J a b2- Ad- 
ditionally, the couplings J c i and J C 2 along the c direction are 
present. 

magnetism pertaining to the pentagonal arrangement of 
magnetic ions in a crystal. 



II. METHODS 

Powder samples of Bi 4 Fe50i3F were prepared by a 
solid-state reaction of Bi 2 3 (99.9 % Aldrich), Fe 2 3 
(nanopowder <50 nm, Aldrich) and BiF 3 (99.99 % 
Aldrich). The initial compounds were weighed in stoi- 
chiometric ratios, mixed, rigorously ground and pressed 
into pellets. The pellets were placed into alumina cru- 
cibles, covered with a lid, and then sealed in quartz tubes 
under a dynamical vacuum of mbar. The samples 

were annealed at 650 °C for 3 h and at 750 °C for 6 h 
with intermediate regrinding. 

X-ray powder diffraction (XRPD) measurements were 
conducted with a Huber G670 Guinier diffractome- 
tcr (CuK a i radiation, curved Ge(lll) monochromator, 
transmission mode, image plate). High-resolution syn- 
chrotron X-ray powder diffraction (SXPD) data were col- 
lected at the ID31 beamline of European Synchrotron Ra- 
diation Facility (ESRF, Grenoble, France) using a wave- 
length of 0.40006 A and eight scintillation detectors, each 
preceded by a Si(lll) analyzer crystal. The powder sam- 
ple was contained in a thin-walled quartz capillary that 
was spun during the experiment. The temperature of the 
sample was controlled by a He-flow cryostat (temperature 
range 10 - 300 K) and hot-air blower (300 - 700 K). 

Neutron powder diffraction (NPD) data were collected 
with the high-resolution powder diffractometer HRPT 
(Paul Scherrer Institut (PSI), Switzerland) at the wave- 
length of 1.8857 A with the use of a standard orange 
cryostat and radiation-type furnace to cover the tem- 
perature range from 1.5 to 700 K. Structure solution 
was performed with direct methods using the EXPO 
software.^ The Rietveld refinement of the crystal and 
magnetic structure against the powder diffraction data 



was performed with the JANA2006 package.^ 

Transmission electron microscopy (TEM) investigation 
was performed with a FEI Tecnai G 2 microscope oper- 
ated at 200 kV. The TEM specimen was prepared by 
grinding the sample in ethanol and depositing the dis- 
persion on a copper grid covered with a holey carbon 
film. The sample was analyzed by means of electron 
diffraction (ED) and high resolution high angle annu- 
lar dark field scanning transmission electron microscopy 
(HAADF-STEM). The simulated HAADF-STEM images 
were calculated using the QSTEM 2.0 softwarePl The Bi 
and Fe content has been measured by energy dispersive 
X-ray (EDX) analysis performed with a Jeol JEM5510 
scanning electron microscope equipped with an INCA 
EDX system (Oxford instruments). The analysis has 
been done using the Bi-M and Fe-K lines. More than 
50 spectra have been collected in order to obtain a sta- 
tistically relevant result. Within the experimental error, 
the measured Bi/Fe = 0.84(6) ratio agrees with the nom- 
inal chemical composition. 

The magnetic susceptibility of Bi4FesOi 3 F was mea- 
sured with a Quantum Design MPMS SQUID magne- 
tometer in the temperature range 2 — 600 K in magnetic 
fields up to 5 T. Heat capacity was measured by a relax- 
ation technique using a Quantum Design PPMS in the 
temperature range 1.8 — 320 K. 

The Mossbauer spectrum was recorded at room tem- 
perature (293 K) in transmission mode using a constant- 
acceleration Mossbauer spectrometer with a nominal 2.02 
GBq Co source in a 6 /xm Rh matrix (conventional 
source). The velocity scale was calibrated relative to a 
30 /zm a-Fe foil using the positions certified for the Na- 
tional Bureau of Standards standard reference ma-terial 
no. 1541; a line width of 0.30 mm/s for the outer lines of 
a-Fe was obtained at room temperature. The spectrum 
collection time was 10 days. The Mossbauer spectrum 
was fitted using the MossA software package.^ 

The electronic structure of Bi 4 Fe50i 3 F was calcu- 
lated within the density functional theory (DFT) frame- 
work using the FPLO code^ and the Perdew-Wang 
parametrization of the exchange-correlation potential™ 
The symmetry-irreducible part of the first Brillouin zone 
was sampled with a k mesh of 72 points for the unit cell 
with 92 atoms. Strong correlations in the Fe 3d shell 
were treated by the DFT+C/ procedure with a Coulomb 
repulsion U = 7 eV and Hund's exchange J = 1 eV.^ 
Magnetic couplings were evaluated by mapping total en- 
ergies of collinear spin configurations onto the Heisenberg 
model. 



III. RESULTS 

A. Crystal structure 

A preliminary XRPD investigation of the Bi 4 Fe50i 3 F 
sample reveals a single-phase material with a tetragonal 
unit cell and the lattice parameters a = 8.30318(5) A, 
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TABLE I. Selected parameters of the Rietveld refinements of the Bi4Fe50i3F structure from the NPD data. 



Formula 
Space group 
T (K) 
a (A) 
6(A) 

z 

Cell volume (A 3 ) 
Calculated density (g/cm 3 ) 
Radiation 
Rf, Rp, RwP 



1.5 
8.28671(9) 
18.0501(2) 

1239.49(2) 
7.188 

0.013, 0.030, 0.040 



Bi 4 Fe 5 0i3F 
P4 2 /mbc 

300 
8.2995(1) 
18.0573(3) 
4 

1243.81(3) 
7.165 

Neutron, A = 1.8857 . 
0.013, 0.029, 0.037 



700 
8.3246(1) 
18.1173(3) 

1255.52(4) 
7.098 

0.018, 0.033, 0.042 



c = 18.0640(1) A. The electron diffraction patterns of 
Bi4FesOi3F are presented in Fig. [2j They can be con- 
sistently indexed on the tetragonal lattice with the unit 
cell parameters as determined from the XRPD data. The 
reflection conditions Okl, k = 2n and hhl, I = 2n suggest 
two possible space groups, P4 2 /m6c and P4 2 bc. The 
structure solution was performed in the centrosymmet- 
ric P^ijmhc space group with direct methods providing 
the positions of all cations and part of the anions in the 
unit cell. The missing anions were located using differ- 
ence Fourier maps. The resulting structural model was 
refined against the SXPD data at the temperatures of 
10 K, 298 K and 673 KP and against the NPD data at 
the temperatures of 1.5 K, 300 K and 700 K (Fig. [3]). 
A magnetic contribution was taken into account for the 
refinement of the crystal structure from the NPD data 
collected at T = 1.5 K (see details in section III-B). Oxy- 
gen and fluorine anions can not be distinguished reliably 
from either X-ray or neutron powder diffraction. There- 
fore, the assignment of the fluorine position at 46 (0,0, 4) 
is based on bond-valence-sum (BVS) arguments. 

BVS arguments strongly favor a complete ordering of 
the O and F anions. The F~ anions linked to four Bi 
cations acquire a BVS of 0.91(2), which is close to their 
nominal valence. For the 2 ~ anions occupying this po- 
sition, the BVS would be as low as 1.20(2), thus demon- 
strating a strong and unrealistic underbonding. The 
anion ordering is additionally supported by the room- 
temperature Mossbauer spectrum that is consistent with 
only two distinct coordination types of Fe 3+ cations (see 
below). 

The refinement of the low-temperature diffraction data 
resulted in low atomic displacements parameters ( ADPs) . 
However, at and above room temperature high ADPs 
of Fl and 01 indicated possible displacements of these 
atoms. Indeed, the displacements of Fl and 01 from 
their special positions reduced the ADPs. This effect is 
mostly of dynamic nature, because at low temperatures 
the displacements are suppressed. The 01 and Fl dis- 
placements observed at high temperatures may be driven 
by subtle adjustments of interatomic distances and coor- 
dination (see BVS analysis in Supplemental Material), 
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FIG. 2. ED patterns of Bi4FesOi3F along main zone axes. 

although further studies of the local structure would be 
necessary to fully elucidate this issue. We have also 
checked the possibility of symmetry lowering. The re- 
finement in the acentric group P4 2 6c neither eliminated 
the disorder, nor reduced the reliability factors. Electron 
diffraction patterns clearly demonstrate the presence of 
glide planes perpendicular to the a and b unit cell axes, 
as well as perpendicular to the diagonals of the ab face. 
Therefore, any further symmetry reduction can be safely 
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FIG. 3. Experimental, calculated and difference NPD pat- 
terns of Bi4FesOi3F after the Rietveld refinement at T = 
1.5 K, 300 K, and 700 K. The bars mark the reflection posi- 
tions. For T = 1.5 K the reflection positions of the magnetic 
supercell are shown. 



chain, thus creating large cavities in the framework. The 
cavities are occupied, by the fluorine-centered tetrahedral 
Bi 4 F groups. 

The structure of Bi4FesOi3F viewed along the c direc- 
tion immediately reveals the pentagonal arrangement of 
the magnetic Fe 3+ cations (compare Fig. [4ji and Fig. [T^,) . 
While in Bi2Fe40g similar pentagonal planes are directly 
stacked on top of each other, Bi4FesOi3F demonstrates 
a more peculiar stacking, with the pentagonal planes in- 
terleaved by the Fe30g octahedra. These octahedra dis- 
rupt the ideal eclipsed configuration within the rutile-like 
chains (Fig. Ep,e). 

The Bi4Fe50i3F crystal structure has been validated 
by HAADF-STEM imaging along the [001] and [100] di- 
rections (Fig. [5|. In HAADF-STEM, the brightness of 
a dot corresponding to the projection of an atomic col- 
umn roughly scales as Z n (1.5 < n < 2), where Z is the 
average atomic number along the column. Therefore, Bi- 
bearing atomic columns appear as the brighter dots and 
form the most prominent contrast in the image. The 
insets in Fig. [5] show HAADF-STEM images calculated 
using the refined atomic positions. The good agreement 
between the experimental and calculated images justifies 
the structure solution. 

The room-temperature Mossbauer spectrum for 
Bi4Fe50i3B221can be decomposed into two paramagnetic 
doublets, contributing to 57(4) % (doublet I) and 43(4) % 
(doublet II) of the total absorption area. Mossbauer pa- 
rameters, such as center shift (CS) and quadrupole split- 
ting (QS), are as follows: 

Doublet I: CS = 0.428(9) mm/s, QS = 0.39(2) mm/s, 

Doublet II: CS = 0.19(1) mm/s, QS = 0.65(2) mm/s. 

The CS values for both doublets evid ence t hat all iron 
in the compound is in the Fe 3+ state! 221241 The CS of 
the doublet I is consistent with Fe 3+ in an octahedral 
coordination™ The significantly smaller CS of the dou- 
blet II is in agreement with the tetrahedrally coordinated 
Fe 3+ (Ref.[24]). The 57(4):43(4) ratio of the octahedrally 
and tetrahedrally coordinated iron corresponds well to 
the anticipated 60:40 ratio in Bi4FesOi3F. 



B. Magnetic structure 



excluded. 

The crystallographic information for Bi 4 FesOi3F is 
summarized in Table [I] Atomic parameters and main in- 
teratomic distances are provided in Tables In] and III 



Further details on the structure refinement are given in 
Supplemental MaterialP2 

The crystal structure of Bi4FesOi3F is illustrated in 
Fig. [4ji-e. It is built of infinite rutile-like chains of edge- 
sharing FeOe octahedra running along the c-axis and 
comprising two Fe positions, Fel and Fe3 (Fig. [4^). The 
Fe2 atoms form pairs of corner-sharing FeC>4 tetrahe- 
dra that link the octahedral chains into a framework. 
These links repeat at every third FeOg octahedron of the 



Below TV = 178 K, extra reflections appear on the 
NPD pattern. These reflections originate from a mag- 
netic ordering, as no changes have been seen in the SXPD 
patterns down to T — 10 K. The magnetic reflections in 
the entire temperature range of T = 1.5 — 178 K can be in- 
dexed with a k = (5,5,0) propagation vector. The mag- 
netic structure at T = 1.5 K was solved using a tetragonal 
magnetic supercell a m = a — b, b m = a + b, c m = c. 
The analysis of possible magnetic symmetries was per- 
formed with the ISODISTORT softwareP The solution 
was found in the magnetic space group Pc^/nW^ 

The magnetic moments of the Fe positions, which are 
symmetrically equivalent in the nuclear structure, were 
kept identical. Although the symmetry does not require 
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TABLE II. Atomic positions 
700 K. 



and displacement parameters (Ui so ) of Bi4Fe50i3F at 1.5 K, 300 K and 



T(K) 



1.5 



300 



700 



Bil, 16i (x, y, z) 

x/a 

y/b 

z/c 

Uiso (A 2 ) 

Fel, 8/ (i,0,z) 

z/c 

u iso (A 2 ) 

Fe2, 8h (x,y,0) 
x/a 

y/b 
u iso (A 2 ) 

01 

Uiso (A 2 ) 

02, 16i (x, y, z) 
x/a 

y/b 

z/c 

U iso (A 2 ) 

03, 8h, x, y, 
x/a 

y/b 
u iso (A 2 ) 

04, 8g (x,x+ l,f) 
x/a 

u iEO (A 2 ) 

05, 16i (x, y, z) 
x/a 

y/b 

z/c 

Uiso (A 2 ) 

Fl 



Uiso (A 2 ) 



0.6693(1) 
0.6594(1) 
0.15718(6) 
0.001 

0.08014(9) 
0.001 

0.8511(2) 
0.8385(2) 
0.001 

4a (0,0,0) 
0.001 

0.2918(2) 
0.8764(2) 
0.5865(1) 
0.001 

0.1347(3) 
0.5907(3) 
0.001 

0.3263(3) 
0.001 

0.5814(2) 
0.8602(2) 
0.3331(1) 
0.001 

8e* (0,0,2) 
z/c = 0.2587(3) 

0.001 



0.6706(2) 
0.6592(2) 
0.15697(7) 
0.0064(3) 

0.0800(1) 
0.0049(4) 

0.8515(2) 
0.8388(2) 
0.0058(4) 

8e* (0,0, z) 
z/c = 0.0072(6) 
0.014(1) 

0.2914(2) 
0.8760(3) 
0.5858(1) 
0.0092(5) 

0.1362(3) 
0.5900(3) 
0.0065(6) 

0.3268(3) 
0.0077(6) 

0.5823(2) 
0.8603(2) 
0.3331(2) 
0.0076(4) 

16i** (a;, 0, z) 
x/a = 0.024(3) 
z/c = 0.2597(7) 
0.007(5) 



0.6713(2) 
0.6592(2) 
0.15676(8) 
0.0195(3) 

0.0798(1) 
0.0173(5) 

0.8516(3) 
0.8390(3) 
0.0201(5) 

8e* (0,0, z) 
z/c = 0.0106(6) 
0.035(2) 

0.2896(3) 
0.8761(3) 
0.5859(2) 
0.0225(7) 

0.1368(4) 
0.5899(4) 
0.0189(8) 

0.3271(3) 
0.0242(9) 

0.5827(3) 
0.8592(3) 
0.3331(2) 
0.0207(5) 

16i" (x, 0, z) 
x/a = 0.022(6) 
z/c = 0.260(1) 
0.04(1) 



Reduced occupancy factors: *g(Fl)— 1 **g(Fl)= | 



*!7(01)= I 



identical magnetic moments for Fel and Fe3, we found 
that these moments are equal within the standard devia- 
tion. Therefore, the corresponding linear constraint was 
introduced without affecting the fit quality. The refined 
magnetic moments at T = 1.5 K are 4.06(6) [Ib for the 
Fel and Fe3 positions and 3.34(6) [Ib for the Fe2 po- 



sition, which is significantly smaller than the expected 
moment of 5 [Ib for Fe 3+ in the high-spin configuration. 
Details on the refined magnetic moment components are 
listed in Table El 

The peculiar non-collinear antiferromagnetic structure 
of Bi4Fe50i3F is illustrated in Fig. [6] All magnetic mo- 
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Fe3 



Fe1 
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Fe1 



FIG. 4. (Color online) The crystal structure of F^FesOisF: clynographic viewing (a); [010] projection (b); [001] projection 
(c); Cairo pentagonal lattice of the Fe atoms (d) and enlarged view on the rutile-like chain (e). FeOe octahedra and Fe04 
tetrahedra are shown in blue and yellow, respectively. The Bi, F, Fe, and O atoms are drawn as large orange, small light green, 
medium green and small blue spheres, respectively. 



merits are confined to the ab plane. The order along 
the rutile-type chains is ferrimagnetic, with two parallel 
magnetic moments on the neighboring Fel atoms, and an 
opposite moment on Fe3. The direction of these moments 
is close to [110] and rotates by 90° between neighboring 
chains, so that the moments in second- neighbor chains 
are opposite and cancel each other. The magnetic mo- 
ments of the Fe2 atoms are ~ 4.2° off the 100 direction 
and form an AFM configuration within each Fe 2 Oy dimer 
(Fig.§. 

The temperature evolution of the magnetic reflections 
is non-monotonic and pinpoints at least two additional 
magnetic transitions that follow the initial AFM ordering 
at T N = 178 K (Fig. [7]). The transition at T x w 62 K is 
manifested by abrupt intensity variations: the reflections 
with even I rapidly vanish, whereas the intensities of the 
reflections with odd I notably increase. The second tran- 
sition at T2 sa 71 K manifests itself by the vanishing of 
the \ I, I — 2n reflections. The two magnetic transitions 
below T/v are well in line with the results of thermody- 
namic measurements (Section IIICl. These transitions 



are likely accompanied by drastic changes in the mag- 



netic structure. Details of this complex and intriguing 
magnetic behavior should be addressed in future studies. 



C. Thermodynamic properties 

The magnetic susceptibility of Bi4Fe50i3F reveals an 
overall antiferromagnetic behavior (Fig. [8]) in agreement 
with the antiferromagnetic, albeit non-collinear spin ar- 
rangement, pinpointed by neutron diffraction. The field 
dependence of the susceptibility should be ascribed to a 
tiny amount of a ferromagnetic impurity, such as Fe2C>3. 
Magnetization loops measured at constant temperature 
do not show any clear signature of a remnant magneti- 
zation or hysteresis (Fig. |9j top). 

Above 300 K, the susceptibility follows the Curie- Weiss 
law x = C/{T + 9) with C = 3.59(3) emu K/mol Fe and 
9 = -380(10) K (Fig. [9) bottom). The resulting Curie 
constant corresponds to an effective magnetic moment 
of /x e ff = 5.4 /is, in reasonable agreement with 5.92 /is 
expected for the spin-5/2 Fe 3+ cations. The negative 9 
indicates predominant AFM interactions in Bi4Fe50i3F. 
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TABLE III. Selected interatomic distances (in A) for 
Bi 4 Fe 5 0i 3 F at 1.5 K, 300 K, and 700 K. 



T(K) 


1.5 


300 


700 


Bil-02 


2.173(2) 


2.173(3) 


2.176(3) 


Bil-04 


2.160(4) 


2.125(2) 


2.133(2) 


Bil-05 


2.170(2) 


2.168(2) 


2.171(3) 


Bil-05 


2.551(2) 


2.555(2) 


2.557(3) 


Bil-Fl 


2.456(4) 


2.35(1) 


2.37(3) 


Bil-Fl 


2.658(4) 


2.57(2) 


2.56(3) 


Bil-Fl 




2.58(1) 


2.61(3) 


rsil— r 1 




2.79(2) 


o on/ o\ 

2.80(3) 






BVS = 3.15(1) 




Fel-02x2 


2.009(2) 


2.017(2) 


2.035(2) 


Fel-03x2 


1.975(2) 


1.981(2) 


1.986(3) 


tel-(J5x 2 


2.064(2) 


2.067(3) 


o no A ( A \ 

2.084(4) 






T)"\ 7"0 O f\/~' A f C> \ 

BVs = 2.964(8) 




Fe2-01 


1.820(1) 


1.824(2) 


1.833(3) 


Fe2-02x2 


1.912(2) 


1.900(2) 


1.900(3) 


Fe2-03 


1.887(3) 


1.882(3) 


1.884(4) 






BVS = 2.93(1) 




Fe3-04x2 


2.035(3) 


2.032(2) 


2.036(3) 


Fe3-05x2 


2.010(2) 


2.015(2) 


2.028(3) 






BVS = 2.956(7) 





TABLE IV. Positions of the magnetic atoms, the components 
of their magnetic moments, and total magnetic moments (M) 
in Bi 4 Fe 5 Oi 3 F at 1.5 K (a = b = 11.7192 A, c = 18.0501 A, 
all m z — 0). 



Atom x/a 



y/b 



z/c 



Fell 
Fel 2 
Fe2i 
Fe2 2 
Fe3 



i 

4 
1 
4 

0.0063 



i 

4 
1 

4 

0.8448 



0.0801 
0.4199 




0.0063 -0.3448 
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In fields of 0.1 T and higher (Fig. |8j bottom part), 
the most conspicuous effect is the double anomaly at 
Ti = 62 K and T 2 = 71 K. In lower fields, a third 
anomaly around T/v = 178 K could be additionally ob- 
served (Fig. M top). Heat-capacity (C p ) measurements 
consistently show the anomalies at Ti, T2, and Tjv, thus 
confirming the intrinsic nature of these magnetic transi- 
tions (Fig. [8j inset). The anomalies at Ti and T2 are in- 
dicative of first-order transitions with small latent heats. 

Considering the neutron diffraction data, we conclude 
that Bi4Fe50i3F undergoes an antiferromagnetic order- 
ing of the Fe moments at T/v — 178 K. However, this or- 
dering is presumably incomplete, as shown by the large 
anomalies in the heat capacity at Ti and T 2 . The area un- 
der the heat-capacity curve plotted as C p /T vs. T corre- 




— Out 

*' ». • • • • m m *MMiilW^l , 'l>4l | i>*'MW * <*» *•» «» 
- • m m m m »» *p mi* * «M»«# •* 

■>* *• ^ m *« «* m» -rm m m «•»■ «• <■* *•* »• ** «* mp 




FIG. 5. [001] and [100] HAADF-STEM images of 

Bi4FesOi3F. The insets show the images calculated using the 
refined crystal structure. The size of the [001] calculated im- 
age is 3a x 3b (thickness t = 4 nm); the size of the [100] 
calculated image is 36 x 2c (thickness t = 7 nm). 



sponds to the entropy release upon a magnetic transition. 
The large amount of entropy released at Ti and T2 im- 
plies that above Ti the system develops an incompletely 
ordered intermediate AFM structure, which is further al- 
tered at T2 and eventually disappears at Tn- This pecu- 
liar behavior is independently confirmed by the drastic 
changes in the magnetic neutron scattering (Sec.III-B). 



D. Microscopic magnetic model 

To elucidate relevant magnetic interactions, the band 
structure of Bi 4 Fe50i3F was calculated. Regarding the 
band structure, this oxyfluoride is very similar toinsu- 
lating Fe 3+ oxides. Its electronic spectrum (Fig. 10 1 re- 
veals oxygen and fluorine 2p bands between —5 eV and 
the Fermi level. The Fe 3d bands are around —6 eV 
(filled) and at +2 eV (empty), with the splitting of 
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FIG. 6. (Color online) The magnetic structure of 

Bi4Fe50i3F at T — 1.5 K. The orientation of the magnetic 
moments within the pentagonal Fel-Fe2 layer at z = \ is 
shown at the bottom. 




50 75 
Temperature (K) 

FIG. 7. (Color online) Temperature variation of the intensity 
of the magnetic reflections. 




100 200 
Temperature (K) 

FIG. 8. (Color online) Magnetic susceptibility of Bi4Fe 5 0i3F 
measured in applied fields of 0.01 T, 0.5 T, and 5 T. The siz- 
able field dependence is due to a tiny amount of ferromagnetic 
impurity. The inset shows the heat capacity divided by tem- 
perature (Cp/T). Arrows denote the transition temperatures 
T\, T 2 , and T N (compare to Fig. 
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FIG. 9. (Color online) Magnetization curve of Bi4FesOi3F 
measured at 2 K (top) and inverse magnetic susceptibility 
(bottom) measured in the applied magnetic field of 1 T (cir- 
cles) along with the Curie- Weiss fit (solid line). 
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FIG. 10. (Color online) Electronic density of states calcu- 
lated for the lowest-energy collinear AFM spin configuration 
of Bi 4 Fe 5 0i 3 F with DFT+U (U = 7 eV). The Fermi level is 
at zero energy. Only one spin direction is shown. 



TABLE V. Magnetic interactions in Bi 4 Fe 5 Oi 3 F. The Ta- 
ble lists Fe-Fe distances, Fe-O-Fe angles, exchange integrals 
Ji, as well as normalized spin-spin correlations {SiSj)/S 2 in 
the ground-state magnetic configuration, as obtained from 
Monte-Carlo simulations. 



Interaction 


d-Fc-Fc (A) 


VFe-O-Fe (deg) 


Ji (K) 


{SiSj)/S 2 


Jcl 


2.91 


94.2 


34 


+0.99 


J C 2 


3.06 


97.4 


10 


-0.94 


Jabl 


3.39 


119.2 


45 


-0.84 


Jab2 


3.53 


130.9 


74 


-0.80 


Jd 


3.64 


180 


191 


-0.97 



8 eV that roughly matches the on-site Coulomb repul- 
sion of U = 7 cV. Bi 6s orbitals form bonding states 
at around —10 eV. In DFT+U calculations, the energy 
gap of 1.5 — 2.0 eV depending on the spin configuration 
concurs with the brown color of Bi 4 Fe50i3F. 

Magnetic couplings extracted from the DFT+U band 
structure indicate a predominant interaction Jd — 252 K 
within the Fe207 dimer, in agreement with the 180° 
Fe-O-Fe exchange pathway, which is most favorable 
for the AFM coupling, according to the Goodenough- 
Kanamori- Anderson rules. Other nearest-neighbor cou- 
plings are also AFM, but much weaker than Jd according 
to the lower Fe-O-Fe angles (see Table [V ) . In general, 
the couplings decrease as the Fe-O-Fe angle approaches 
90°, although this trend is violated at the lowest angles 
(J c i > J C 2)i where the direct Fe-Fe exchange comes into 
play. The couplings beyond nearest neighbors are all be- 
low 2 K and can be safely neglected in a first approxima- 
tion. The resulting spin lattice has pentagonal-like units 
and 3 nonequivalent couplings in the ab plane (Fig.JTJ)). 
The interactions J C 2 connect these planes along the c di- 
rection via the Fe3 spins, which do not belong to the 
pentagonal units. 

Monte-Carlo simulations for the spin lattice depicted in 



Fig. [TJa and the J, parameters from Table |V| suggest that 
at low temperatures Bi4Fe50i3F should be magnetically 
ordered. The ordering pattern can be extracted from the 
spin-spin correlations (Table |v|) , as explained in Ref. 121)1 
The normalized correlation (SlSj)/S 2 is the measure of 
cos (p, where tp is the angle between the magnetic mo- 
ments. Negative normalized correlations approaching —1 
indicate antiparallel magnetic moments within the Fe20y 
dimers (the interaction Jd) as well as between Fel and 
Fe3 along the c direction (the interaction J C 2)- Despite 
the sizable AFM coupling J cl between the adjacent Fel 
sites along c, the moments on these Fel sites are parallel, 
as shown by the positive spin-spin correlation of nearly 
+1 (Table [v]). In this way, the simulations identify the 
ferrimagnetic order within the rutile-like chains, in agree- 
ment with the experiment. 

Following the expectations, the spin lattice of 
Bi 4 Fe50i3F features pentagonal units with exclusively 
AFM interactions (Fig. [T]d) that lead to a strong mag- 
netic frustration. Indeed, in the ground-state magnetic 
structure only the strongest coupling Jd is fully satisfied. 
The moments coupled by J a bi and J a hi develop the non- 
collincar configuration, whereas the AFM coupling J cl 
is fully suppressed and gives way to the FM order. The 
coupling J C 2 is also fully satisfied, because it does not 
compete with any other interaction in this system. 



IV. DISCUSSION AND SUMMARY 

Our experimental and theoretical results establish the 
strong magnetic frustration in Bi 4 Fe50i3F. This effect 
originates from the pentagonal units of the spin lattice. 
The frustration underlies the intricate non-collinear mag- 
netic structure and the overall peculiar magnetic behav- 
ior with the magnetic ordering at Tjv and further trans- 
formations at T\ and T 2 . 

Similar to Bi2Fe 4 0g, Bi 4 FesOi3F is not a perfect pro- 
totype of the Cairo spin lattice (Fig. [I]). In Bi 4 FesOi3F, 
two nonequivalent couplings J a bi and J a bi stand for the 
single coupling Ji- Additionally, the lattice is decorated 
by Fel dumbbells that replace all four-vertex sites. How- 
ever, this kind of decoration does not change any ba- 
sic features of the magnetic system. The spin lattices of 
Bi2Fe 4 0g and Bi 4 FesOi3F are topologically equivalent to 
the Cairo lattice. The non-collinear magnetic structures 
observed in these compounds follow theoretical predic- 
tions for the ideal Cairo lattice.^ 

Two striking differences between Bi 2 Fe 4 0g and 
Bi 4 FesOi3F should be emphasized, though. Experimen- 
tally, Bi2Fe 4 09 orders antiferromagnetically below 264 K 
and retains the same magnetic structure down to the low- 
est temperatures. This can be seen from the s mooth 
temperature evolution of the magnetic reflections! 14 * 27 ! 
In Bi 4 Fe50i3F, the AFM ordering that sets in below 
Tjv = 178 K is likely incomplete and changes upon fur- 
ther magnetic transitions at T\ and T 2 . As none of 
these transitions is accompanied by any structural trans- 
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formations detectable by NPD and SXPD, the sizable 
release of entropy and the indications for a weak first- 
order contribution imply that above T\ the magnetic or- 
der in Bi4FesOi3F is altered. One plausible explana- 
tion is the formation of partially disordered AFM states 
that have higher entropy than the fully ordered non- 
collinear AFM structure. The partially disordered states 
are found in magnetic models based on triangular spin 
latticesP^Hl They are indeed observed in Co 2+ halides, 
such as CsCoCbJ^and RbCoBr 3 . 32 Surprisingly, our ex- 
perimental results may put forward similar physics on the 
pentagonal lattice. Our results also rise further questions 
to theory, because the sequence of magnetic transitions 
below the Neel temperature was never anticipated for any 
pentagonal geometry. 

The second difference between P^FesOisF and 
Bi 2 Fe40g pertains to the stacking of the pentagonal lay- 
ers. While in Bi 2 Fe40g these layers are directly stacked 
on top of each other, Bi4FesOi3F features additional Fe3 
sites between the pentagonal layers. The formation of 
these additional sites accompanies the incorporation of 
the Bi 4 F tetrahedra. Indeed, a comparison between the 
crystal structures of Bi 2 Fe40g, Bi 4 FesOi3F, and silli- 
manite Al 2 SiOs (the aristotype structure for the mullitc- 
type compounds^) opens broad perspectives for finding 
new systems with a pentagonal arrangement of magnetic 
atoms. All these structures contain infinite rutile-type 
chains of octahedra and differ by the connectivity of 
the chains through the B2O7 tetrahedral groups. In the 
Al 2 Si05 (= B3O5) sillimanite structure, the connectiv- 
ity is the densest: every BO@ octahedron shares a corner 
with the tetrahedral group (Fig. 11). In the Bi 2 Fe40g 
structure, the linkages are repeated at every second oc- 
tahedron, whereas in Bi 4 Fe 5 0i3F every third octahedron 
is linked. 

The series of the Al 2 Si0 5 , Bi 2 Fe 4 9 , and Bi 4 Fe 5 0i 3 F 
compounds can be formally considered as an insertion 
of A 2 B04 fragments between the tetrahedral linkages. 
The A 2 B04 entity includes two atomic layers, A 2 2 and 
B0 2 , and increases the distance between the tetrahe- 
dral linkages. Additionally, two A-cations located be- 
tween the chains appear. This way, one arrives at the 
sillimanite-based homologous series B3O5 + nA 2 B04, 
where n — corresponds to the sillimanite structure, 
n = 1 compri ses a gr oup of Bi 2 Fe 4 O g (B = Fe, Al, 
Ga) compound j 12 l 13 l 3 4ri3H ancl n — 2 is represented by 
Bi4Fe50i3(F), where the extra F atom is added between 
two A 2 2 layers. Extending this series ton = 00, one ar- 
rives at the A 2 B04 structure, where tetrahedral linkages 
between the rutile-type chains are absent. This structure 
is exemplified by Sb 2 B0 4 (B = Mn, Fe, Ni^Efll belong- 
ing to the Pbs04 structure type. 

The above crystallographic arguments suggest that 
other members of the homologous series should be formed 
in between the limiting cases of B3O5 (n = 0) and A 2 B04 
(n = 00). Indeed, the n = 1 (Bi 2 Fe40g) and n = 2 
(Bi4FesOi3F) members are already available. The in- 
termediate compounds would contain pentagonal-lattice 



n = 1 




n = 2 



OOOOO n=oo 




FIG. 11. (Color online) The scheme of the linkage between 
rutile-type octahedral chains (blue) through the tetrahedral 
units (yellow) in the B3O5 + nA2B04 homologous series for 
different n. The A cations are shown as orange spheres. 



units in the ab plane connected via intermediate B sites 
along the c direction. While Bi 2 Fe40g shows a single 
magnetic transition, Bi4FesOi3F reveals a more intricate 
and peculiar physics with several transitions that can be 
retained or probably altered in further members of the 
family with higher n. This approach also enables the 
control over the dimensionality, because the separation 
between the pentagonal layers increases with n. 

The compounds with high n require appropriate 
A cations. Considering the available Bi 2 Fe 4 0g and 
Bi4FesOi3F compositions, the lone electron pair is likely 
a necessary prerequisite. Besides Bi 3+ , Sb 3+ and Pb 2+ 
are plausible candidates because of their known com- 
patibility with the A 2 B04 structurePS The combina- 
tion of Bi 3+ , Sb 3+ and Pb 2+ enables heterovalent re- 
placements and the fine tuning of the electron concen- 
tration. Another interesting feature is the ability to 
accommodate extra anions. The F _ anions can be in- 
serted into tetrahedral voids formed between two succes- 
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sive A2O2 layers, as in the n = 2 Bi4Fe50i3F structure 
that has the A4B5O13F composition instead of the ex- 
pected A4B5O13. Alternatively, the anion content can 
be changed by replacing the common vertex of two Fe04 
tctrahedra with an oxygen dumbbell that would trans- 
form a pair of corner-sharing tetrahcdra into a pair of 
edge-sharing tetragonal pyramids B 2 Os- This replace- 
ment is observed in the Mn-based version of the n = 1 
compound, Bi 2 Mn 4 Oi (= A 2 B 4 O g O)!±2] 

In summary, the discovery of Bi4Fe50i3F not only pro- 
vides an interesting material prototype of a frustrated 
Cairo pentagonal lattice, but also opens tantalizing op- 
portunities for a further expansion of this interesting 
class of materials. At low temperatures, Bi 4 FesOi3F de- 
velops the non-collinear antiferromagnetic order, which 
is anticipated for the pentagonal Cairo lattice. How- 
ever, this order is strongly altered at T\ sa 62 K and 
T 2 ~ 71 K. The respective magnetic transitions have no 
detectable structural component. They also involve the 
sizable increase in the magnetic entropy that may be as- 
cribed to the transformation toward partially ordered an- 
tiferromagnetic states extending up to the Neel tempera- 



ture of Tjv ~ 178 K. Additionally, Bi4Fe50i3F manifests 
the B3O5 + nA 2 B04 homologous series of compounds 
that form further examples of pentagonal spin lattices, 
thus facilitating a precise tuning of the "dimensionality" 
as well as electronic concentration in these systems. 
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